Vortices and quasiparticles near the "superconductor-insulator" transition in thin films 
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We study the low temperature behavior of an amorphous superconducting film driven normal by a perpendicu- 
lar magnetic field (B). For this purpose we introduce a new two-fluid formulation consisting of fermionized field 
induced vortices and electrically neutralized Bogoliubov quasiparticles (spinons) interacting via a long-ranged 
statistical interaction. This approach allows us to access a novel non-Fermi liquid phase which naturally interpo- 
lates between the low B superconductor and the high B normal metal. We discuss the transport, thermodynamic, 
and tunneling properties of the resulting "vortex metal" phase. 



o 
o _ 

u • 
Of 

3 ; 



i 

-a 
c 

o 
o 



(N 
> 

in 

o 
in 
o 

-I— > 



i 

C 

o 
o 



X 



Superconductivity in two dimensions (2d) provides a 
unique area to explore a fascinating variety of quantum phe- 
nomena. Of particular interest are strongly disordered super- 
conducting films, tuned by a perpendicular magnetic field into 
the normal state. Early experiments | lj yl focused on a 
magnetic-field tuned "superconductor-insulator" transition at 
a critical field B c . Based on a phenomenological "dirty bo- 
son" model |4], universal scaling behavior in the temperature 
and field dependence of the resistance was predicted near the 
quantum phase transition (B « B c ) separating the supercon- 
ductor from the Bose insulator. Although some experimental 
support was found for scaling of R(T, B), the resistance at the 
separatrix, R(T - 0, B c ) was often found to deviate signifi- 
cantly from the expected universal value (near Rq = h/4e 2 ). 

More recent experiments on such amorphous films at lower 
temperatures revealed a rather rich and surprising behavior in 
R(T — > 0, B). Specifically, for temperatures well below the 
mean field transition where Cooper pairs are well formed, the 
resistance saturates at a value which can be many orders of 
magnitude smaller than the normal state resistance |5, 6J. This 
mysterious strange metal phase occurs over a reasonably large 
range of intermediate fields, B < B c . Moreover, in highly dis- 
ordered films with B c <K H C 2, the resistance at higher fields 
was found to exhibit a dramatic non-monotonic dependence, 
rising in some cases up to values of order 10 6 /?g, before plum- 
meting towards the normal state resistance at B « H C 2 LZlLa-bS] • 
A boson-only theory that might exhibit a metallic phase was 
proposed by Dalidovic et al. fiolfllll . Kapitulnik et al. fi^l . 
however, argued for the importance of fermionic quasiparti- 
cles ignored in the dirty-boson approach. 

In this Letter we develop a two-fluid formulation, which 
incorporates fermionic quasiparticles into the dirty-boson 
model. The two fluids consist of field induced vortices and 
electrically neutralized quasiparticles ("spinons") interacting 
via a long-ranged statistical interaction. In the low temper- 
ature limit the vortices must be treated as quantum particles 
and their statistics will be important. Previous work has im- 
plicitly assumed that the vortices are bosonic. However in 
the presence of unpaired electronic excitations the issue of 
vortex statistics can be subtle. In this paper we explore the 
consequences of treating the vortices as fermions (see also 
Il3l) . Technically this is achieved via a Chern-Simons flux at- 
tachment to the bosonic vortices. With fermionic vortices and 



fermionic electron-like quasiparticles our two-fluid model can 
be studied within an RPA approximation. The fermionized 
vortices can quantum diffuse, leading to a description of a 
novel conducting non-Fermi liquid phase which interpolates 
between the low field superconductor and the high field nor- 
mal state. We suggest that this "vortex metal" phase might 
account for the strange metallic behavior observed in InO 
and MoGe films. Moreover, our two-fluid approach gives a 
natural explanation for the non-monotonic dependence of the 
low temperature magnetoresistance. Below we discuss trans- 
port, thermodynamic, and tunneling properties of the "vortex 
metal." 

Consider a 2d disordered superconductor in a perpendic- 
ular magnetic field below the mean-field BCS upper critical 
field H C 2- For concreteness we consider a lattice tight-binding 
Hamiltonian of an s-wave superconductor, 



( H=-t^Y_ k [c\ a c? a e iK "' + h. c] + — Yj ( 2 "r + c\ a c Ia - n ) 
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where A is the magnitude of the superconducting order param- 
eter, C~ l is an on-site charging energy and no is the electrical 
charge density. At temperatures well below the mean-field 
transition, it is necessary to include quantum fluctuations of 
the Cooper pairs, which were incorporated phenomenologi- 
cally via the phase of the superconducting order parameter, 
4>i, which is conjugate to the Cooper pair number operator, n r . 
Due to the external magnetic field and strong disorder, one ex- 
pects a substantial number of low energy electrons will also be 
present, certainly in the vortex cores and perhaps elsewhere. 
Moreover, since the strange metal is resistive, the vortices are 
mobile, and a correct description will likely require incorpo- 
rating the quantum dynamics of both vortices and fermionic 
quasiparticles. Such a formulation has recently been obtained 
in Refs. fl4l[l5ll . and will be reviewed and exploited herein. 
Following Ref. 1 14] it is convenient to make a change of 

variables: f} a = b T c\ a , where (t>l^j = e'^ r creates a Cooper 
pair at site r. This transformation removes the electric charge 
from the quasiparticles, leaving charge zero and spin one-half 
fermionic "spinons", f r . The use of spinon variables is merely 
a technical convenience and does not necessarily imply that 
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they are the good excitations in the system. Indeed in the vor- 
tex metal phase the spin-carrying excitations will ultimately 
be electrons (not spinons). Next, a duality transformation can 
be implemented which exchanges the operators b r for hc/2e 
vortices, leading to a theory of bosonic vortices in a fluctu- 
ating gauge field which mediates the inter-vortex long-range 
interactions. Moreover, the vortices and spinons have a statis- 
tical interaction between them, i. e., a vortex "sees" a spinon 
as a source of a 7r-flux and vice versa. 

We describe the vortex-spinon mixture with a so-called 
U(l) formulation fl5l Hill . In the long-wavelength low- 
energy limit the corresponding Euclidean action is: S - 
J dxdT[X, v + X s + -CintL where X v is the vortex Lagrangian 
in terms of the vortex bosonic fields *F and the gauge field 



(a, fl ) = (a o a lt , a ), 



Xv = 



( - ia + ia) 

2m v 
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¥ + Za, (2) 



with: Xa — j£ (e^dyd^ - d^nno} . The spinon's action is 



-£s — f a 



2m, 



(V - i/3) z + (d T - iB ) 



fa, 



(3) 



and the vortex-spinon statistical interaction is mediated by two 
auxiliary f/(l) fields andyS^, 



X, 



-{il^e^Aa^dvPA 



(4) 



The equations of motion S£/Sao = and 6£,/6Bq = attach 
ji flux to the vortices and spinons, 

etjdfij = nN v = n^Sf; e^aj = nh { = nfj a . (5) 

with N w and hf the vortex and spinon densities. The total elec- 
trical charge density is given by no = (ey<9,a/)/7r, and for small 
capacitance will be set by the c-number no, that is (no) = »o- 
The average number of vortices is set by the external magnetic 
field B through n (N v ) = B (in units where fi — c — e — 1). 

In the strange metal, where the film resistance R(T) satu- 
rates at low temperatures, it appears that the vortices are dif- 
fusing. To access this within the vortex-spinon theory, it will 
prove extremely convenient to statistically transmute the vor- 
tices, converting them into fermions. This can be achieved 
by attaching 2n "statistical" flux to each vortex, introducing a 
Chern-Simons gauge field A^. If we denote the Lagrangian in 
Eq. Q as X v ( v f, cifi, the Lagrangian for the fermionized 
vortices (denoted as if/) is simply, 



Xf v = a-A,a)+ —A^e^dyA^. 



(6) 



As recently argued, fermionization of vortices is expected to 
be an extremely good approximation due to the long-ranged 
intervortex interaction 1 1 711 - Indeed, by defining a'^-a^-A^, 
and then integrating over Ap, we absorb the Chern-Simons 
gauge field A^ into a^. One thereby obtains, 



Xfv =£ v ty,a',a) + 0(d 3 a' 2 ) 



(7) 



where we have dropped terms less important than the Maxwell 
term present in £ y . Remarkably, the resulting theory of 



fermionized vortices has no Chern-Simons term. The full La- 
grangian, + X s + Xint> describes fermionic vortices and 
spinons interacting via a statistical interaction, and constitutes 
our two-fluid formulation of 2d superconductors in a field. 

Generally, the fermionized vortices will see an effective av- 
erage "dual" magnetic field with strength 



bfy - n(ho - hf - 2N V ). 



(8) 



The electrical charge density and the vortex density are exter- 
nally determined conserved quantities. But now consider hf. 
as it stands, in addition to the conservation of total electrical 
charge, the total spinon (or electron) number, Nf = J dx(hf) 
is also conserved by X- This additional global symmetry is 
present because we have dropped a term proportional to A, 



Xa = A[/ rT / ri (?!(r) + h.c.l 



(9) 



where the operator (^(r) inserts 2n flux in eijdjaj at r. This 
term originates microscopically from the term e^'^cc whereby 
a Cooper pair is created from two electrons. As we will later 
show, the term Xa is perturbatively irrelevant in the "vortex 
metal" phase. Thus the conservation of Nf is an "emergent 
symmetry" not present microscopically. Although the term 
in Eq. © is irrelevant, it has an important role: it allows the 
spinon density to adjust, such that the total energy of the sys- 
tem with A = is minimized. This condition then determines 
the spinon density hf. It is possible that over some range of 
parameters the spinon density adjusts itself to make bf v = 
and fermionic vortices are natural variables. Another example 
is the limit of vanishing vortex mass, m v — > 0. In this case, 
the bosonic vortices "condense" and expel the flux (b v = 0), 
which sets (hf ) = no, effectively "gluing" the charge back on 
to the spinons. Thus, one recovers an ordinary Fermi liquid 
of electrons. Alternately this may be described as an integer 
quantum Hall state for fermionized vortices with v v = 1 . 

Despite the average effective field in Eq. 0, the transverse 
force on a moving vortex is expected to be small in the dis- 
ordered situations of interest - essentially due to the "normal" 
core, the local dual magnetic field in the vicinity of the vortex 
is small. In this case, we expect that the fermionized vortices 
can quantum diffuse, leading to a non-zero resistance at low 
temperatures. Below we explore the properties of the "vor- 
tex metal" phase resulting from the action in Eqs. J2I6I >. We 
neglect possible quantum interference effects. 

For transport properties, a semi-classical Boltzmann equa- 
tion or even Drude theory is probably adequate. In the lat- 
ter case, the statistics of the particles are not important, but 
it is essential that the statistical interaction between vortices 
and spinons is correctly incorporated. Specifically, when the 
spinons move they induce an electro-motive force (EMF) on 
the vortices and vice versa. Similarly, an electrical charge cur- 
rent induces an EMF on the vortices. These effects can be 
summarized in the following transport equations, which de- 
scribe the response of the system to an electrical current J and 
to a thermal gradient WT: 



f v = o- v -zx(j s - J)-i v vr ; 

J s = cr s • z X J v - A S VT, 



(10) 
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FIG. 1 : Schematic T = phase diagram of amorphous superconduct- 
ing films in the two-fluid model, as a function of an applied magnetic 
field. Dashed curves label the vortex conductivity and spinon resis- 
tivity, which determine via Eq. 1 1 1 i the electrical resistance (solid 
line). The "vortex metal" phase separates the low field superconduc- 
tor (B < B 2 ) from the high field Fermi liquid (B > B 3 ). 



where j v / s are vortex and spinon currents, respectively, and z 
is a unit vector normal to the film. Here, <x v / s and /t v / s are 
dimensionless transport "conductivities" of the bosonic vor- 
tices/spinons determining their current response to the EMF 
and thermal gradient. Generally these quantities are two-by- 
two matrices with off-diagonal components. The resistivity 
tensor of the fermionic vortices (<Xf v ) will differ from those 
of the bosonic ones by integer off-diagonal terms due to the 
Chern-Simons flux attachment. Thus, (cr fv ) _I = (cr v ) _I + 2ne, 
with e the unit antisymmetric tensor (e xy = 1). The physical 
electrochemical potential gradient is given by £ — z x j v . Note 
that both the magnetic field and the electric field seen by the 
spinons is the same as the external ones. Thus the spinons re- 
spond to the external electromagnetic field as electrons. In the 
absence of a thermal gradient one can extract the (dimension- 
less) electrical sheet conductance matrix defined by, J = cr&. 



(cr v ) 1 + cr s 



(ID 



Thus the spinon and vortex contributions to the total conduc- 
tivity simply add within this two-fluid theory. It is expected 
that the off-diagonal components of cr~ 1 and <x s are both small 
compared to the diagonal component. Ignoring them we ob- 
tain for the longitudinal resistance R = o~$ x /(l + o- xx cr xx ). 

Before exploring the "vortex metal" phase, we discuss the 
expected (T = 0) behavior of <x v / s as the strength of the exter- 
nal magnetic field is varied (Fig.[Q. At low field the vortices 
are well separated and localized by impurities and the strong 
inter-vortex interaction, implying cr xx = 0. This is the super- 
conducting (vortex glass) phase with R(T — 0) = 0. With 
increasing field and density, the vortices can delocalize at a 
field denoted by B2 in Fig. [2 and will quantum diffuse with 
nonzero diagonal cr xx and nonzero electrical resistance - this 
is the "vortex metal." Finally, at large fields (B3) the bosonic 
vortices can condense giving cr xx = 00 (the fermionized vor- 
tices form the v = 1 IQHE state with cr xx = 0). This is the 



conventional Fermi liquid with R - l/<x s . The spinon con- 
ductivity is likewise expected to be zero at very low magnetic 
fields, since the spinons will predominantly be localized at 
the cores of the well separated vortices. At some field, B\, 
spinons may become delocalized, and form an impurity band 
connecting vortices |18]. If this occurs with cr xx = 0, this is 
a "superconducting thermal metal" phase. Finally, in the very 
high field limit, one expects the spinon resistivity to approach 
the normal Fermi liquid resistance, (cr xx )~ l = R N . 

For very disordered films with weakened superconductiv- 
ity, the vortices will be especially mobile, and will condense 
at low magnetic fields, B3 <k H C 2- With dilute vortices 
the spinon resistivity could be significantly larger than R^ at 
B « B3, and the films electrical resistance will have a very 
large peak just below B3. This behavior is consistent to that 
observed in such disordered samples. 

Let us now consider the thermal conductivity k xx . We as- 
sume that Fermi-liquid theory applies separately to both the 
fermionized vortices and spinons, then k xx = T(cr xx + cr xx ) 
(setting the Lorentz ratio to be 1 for each fluid). Thus the 
Lorentz ratio for the two-fluid vortex metal is L = k xx /<t xx T = 
(cr xx + o- xx )/[<t xx + (<t xx )~ 1 ]. This violation of the Wiedemann- 
Franz law is a striking signature of the non-Fermi liquid na- 
ture of the vortex metal. Notice that L varies from zero to 
one as the magnetic field is tuned through the "vortex metal" 
phase from B2 up to B3. Thus, the "vortex metal" naturally 
interpolates between the superconducting and normal phases. 
As for the Nernst effect, N = & y /d x T, with 7 = 0, ignoring 
small contributions from the off-diagonal conductivities we 
find, N = {A xx l<j xx )R. Within Fermi liquid theory (applied to 
the fermionic vortex fluid) A xx oc Tcr xx , the Nernst signal is 
predicted to exhibit a non-monotonic B-field dependence. 

To address the Hall response, we need to consider non-zero 
off-diagonal "conductivities". The Hall angle is then, 



tanf9 H 



ov + Act, 



A = (O 2 + (of) 2 - (12) 



cr xx + Act" 

Another quantity of interest is the Nernst angle. Again ignor- 
ing small off-diagonal contributions, we get 



tan#N 



(13) 



Thus the transport angles 0u, #n differ in the "vortex metal". 

Next we consider the specific heat, C(T), in the "vortex 
metal" at low temperatures. Naively, the vortex and spinon 
Fermi seas would give the usual C ~ T behavior. But this 
is modified due to gauge fluctuations. Within RPA we inte- 
grate out the (diffusive) fermionic spinons and vortices to ob- 
tain the quadratic action, S e ff = S a + Sm + Srpa with, Srpa = 
j qajn \p v \a -a Q \ 2 + D v \a t -a,\ 2 +p s \B Q \ 2 + D,\B r \ 2 ], where a t , a t , 
and B t denote transverse gauge fields in the Coulomb gauge, 
p v /s are the vortex/spinon density of states and £) v / s (q, w n ) = 
0"v/s|w„| +Xv/sq 2 , with "magnetic susceptibilities" Xv/s- The 
free energy can now be readily extracted from the effective 
quadratic action giving, C(T) ~ T In {I IT) + 0(T). 

Let us now discuss tunneling properties. Electron tunneling 
from a metallic tip is expected to be similar to that in a con- 
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ventional 2d metal due to the gapless electron-like quasiparti- 
cles present in the vortex metal state. With a superconducting 
tip placed above the film, one can measure the Cooper pair 
tunneling conductance into the vortex metal 1 19]. The con- 
ductance G = Ga + G cp has two parts: Ga is due to simultane- 
ous tunneling of two electrons, which is equivalent to the An- 
dreev reflection of an incident hole, and G cp is due to tunnel- 
ing of Cooper-pairs that get absorbed in the Cooper-pair fluid. 
Whereas Ga measures the local density of states of the gap- 
less quasiparticles and is almost temperature independent, G cp 
depends strongly on temperature. G cp is given by the Kubo 
formula, G cp oc i J dtt([I(t), 1(0)]). With the tip at x = 0, 
the Cooper-pair current operator is / = 2e/sin0o- First, we 
evaluate the imaginary time correlator, C(t) = /e^Me - ^ )^. 
The density of Cooper pairs is the "dual magnetic field," 
b - eijdiOj and the Cooper pair current is z x e with "electric 
field", e = — V«o - <5oa. C(r) describes the process of insert- 
ing a monopole at x = and removing it at a later time t. 
The action of this monopole configuration determines C(t); 
we calculate it within a quadratic action for the gauge field. 
Possibly, higher order terms in the gauge field action affect 
the result, but they are beyond the scope of this paper. Since 
the spinons do not affect the dominant large r behavior 1 20], 
we focus on the action due to vortices. Integrating them out 
within RPA, we obtain 



T v - " 1 C r 21 
S ° = 7 Zj 7w [B a pe a e p + I ib \ + S /t 



(14) 



So is expressed using the gauge invariant fields e and b. 
Here S ^ — i J dxA(x) [V x e + dob - p m (x)] enforces the dual 
Maxwell equation, with A(x) a Lagrange multiplier and p m (x) 
the density of monopoles; we set p m (xo, r) = 5(xo)d(r)—6(xQ— 
T)6(r). In Eq. (I14> . s a p is the dielectric constant determined 
by the diffusive fermionic vortices, with transverse part = 
cr v /|w„|. The permeability jj — X\ is the vortex orbital "mag- 
netic" susceptibility. Using Eq. fl!4i we calculate the action 
of the monopole configuration and find, S (t) = In 2 (t/to), 
with r ' the scattering rate of the diffusive vortices. Inserting 
the correlator C(j) oc exp [— S(t)] into the Kubo formula and 
performing analytical continuation we find: 



'cp 



(T) oc T 1 exp 



-^ln 2 
2 



(Tt ) 



(15) 



Note that the Andreev part of the conductance, Ga, which de- 
scribes two electron tunneling, adds a ^-independent contri- 
bution to G, which traces <x s : Ga 00 o" 2 . Near the su- 
perconducting transition, Z?2, G cp is expected to dominate the 
tunneling since cr v is small. Particularly, when vortices are lo- 
calized G cp diverges, which corresponds to the appearance of 
the Josephson effect. As the magnetic field, cr v , and cr s grow, 
Ga will increase and the temperature dependence of G gets 
suppressed. Therefore, tunneling provides another probe of 
the superconductor-insulator cross-over in the "vortex metal." 

We can finally address the effects of Xa in Eq. (|9} in the 
"vortex metal" within the RPA. The perturbative effects of Xa 
can be extracted from the two-point correlation function of 



f\f\O a evaluated at A = 0. Equivalently we can consider 
the correlators of e'^CjCj. Assuming factorization to an elec- 
tron correlator times a Cooper pair correlator, the electrons 
contribute a power law since (c a (r)c£(0)) ~ t _1 . The imagi- 
nary time correlator is C(t) oc exp f— cr v ln 2 (r/ro)J. Since this 
decays faster than any power, we expect that in the "vortex 
metal" a small A is perturbatively irrelevant in the RG sense. 

In this Letter we introduced a new two fluid description 
for amorphous superconducting films in a magnetic field. 
Fermionization of the field-induced vortices allowed us to ac- 
cess a novel resistive quantum-liquid phase that naturally in- 
terpolates between the low B superconductor and the normal 
metal for B > H C 2- Anomalous thermoelectric transport and 
tunneling behavior were predicted in the "vortex metal." It is 
our hope that this paper will help motivate new experiments 
on such 2d amorphous films. 

It is a pleasure to thank A. Auerbach, L. Balents, A. Kapit- 
ulnik, S. Kivelson, A. I. Larkin, O. Motrunich, Y. Oreg, and 
D. Shahar for discussions on this work. This work was gener- 
ously supported by the David and Lucile Packard foundation 
(V.G.) and by the National Science Foundation through grants 
PHY99-07949 (V.G., G.R., and M.PA.F), DMR-0210790 
(M.P.A.F), and DMR-0308945 (T.S.). T.S. also acknowl- 
edges funding from the NEC Corporation, the Alfred P. Sloan 
Foundation, and an award from The Research Corporation. 



[1] D. B. Haviland, Y. Liu, and A. M. Goldman, Phys. Rev. Lett. 
62, 2180 (1989). 

[2] A. F. Hebard and M. A. Paalanen, Phys. Rev. Lett. 65, 927 
(1990). 

[3] V. F. Gantmakher and et al., JETP Lett. 68, 363 (1998). 
[4] M. P. A. Fisher, Phys. Rev. Lett. 65, 923 (1990). 
[5] N. Mason and A. Kapitulnik, Phys. Rev. Lett. 82, 5341 (1999). 
[6] N. Mason and A. Kapitulnik, Phys. Rev. B 64, 060504 (2001). 
[7] M. A. Palaanen, A. F. Hebard, and R. R. Ruel, Phys. Rev. Lett. 
69, 1609 (1992). 

[8] G. Sambandamurthy, L. W. Engel, A. Johansson, and D. Sha- 
har, Phys. Rev. Lett. 92, 107005 (2004). 
[9] M. A. Steiner, G. Boebinger, and A. Kapitulnik, Phys. Rev. Lett. 
94, 107008 (2005). 
[10] D. Dalidovic and P. Phillips, Phys. Rev. B 64, 052507 (2001). 
[11] D. Dalidovic and P. Phillips, Phys. Rev. Lett. 89, 027001 
(2002). 

[12] A. Kapitulnik, N. Mason, S. A. Kivelson, and S. Chakravarty, 

Phys. Rev. B 63, 125322 (2001). 
[13] M. V. Feigelman, V. B. Geshkenbein, L. B. Ioffe, and A. I. 

Larkin, Phys. Rev. B 48, 16641 (1993). 
[14] T. Senthil and M. P. A. Fisher, Phys. Rev. B 62, 7850 (2000). 
[15] L. Balents and M. P. A. Fisher, Phys. Rev. B 71, 0851 19 (2005). 
[16] N. Read, Phys. Rev. B 58, 16262 (1998). 
[17] J. Alicea, O. I. Motrunich, M. Hermele, and M. P. A. Fisher, 

cond-mat/0503399. 
[18] S. Vishveshwara, T. Senthil, and M. P. A. Fisher, PRB 61, 6966 

(2000). 

[19] O. Naaman, W. Teizer, and R. Dynes, Phys. Rev. Lett. 87, 
097004 (2001). 

[20] G. Refael, V. M. Galitski, and M. P. A. Fisher, in preparation. 



